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I. INTRODUCTION 



The inflationary paradigm has been remarkably successful in predicting observed fea- 
tures of the cosmic microwave background (CMB) radiation and the large scale structure 
of the Universe. If inflation is driven by a nearly free, massless quantum field, then a 
generic prediction is a spectrum of primordial fluctuations which is Gaussian and almost 
scale- invariant l|-|5| . For a recent review, see for example 6| . The best test of these predic- 
tions comes from CMB observations by the WMAP satellite, which has found a spectrum 
of temperature fluctuations consistent with Gaussian, nearly scale-invariant primordial fluc- 
tuations [7|. 

However, in addition to the dominant effect coming from tree order fluctuations of the 
scalar inflaton and of the graviton, there should also be some effects from loop corrections of 
these fields with themselves and with other fields. The latter will be the topic of this paper, 
particularly a one loop effect which can be interpreted in terms of quantum stress tensor 
fluctuations. The fluctuations of quantum stress tensors and their physical effects have 



been discussed by several authors in recent years 



12 1 . For a recent review with further 



references, see Ref. 13|. Quantum stress tensor fluctuations necessarily have a skewed, highly 



non-Gaussian, probability distribution, although the explicit form of this distribution has 



only been 
Ref. 



bund in two-dimensional spacetime models 14j |. 
15| studied the possible contributions of quantum stress tensor fluctuations of a 
conformally invariant field to primordial density perturbations in inflationary models. It 
was found that these contributions can be proportional to a power of the scale factor change 
during inflation, and hence potentially large enough to observe. Because they are asso- 
ciated with a non-scale invariant and non-Gaussian contribution, they can at best be a 



sub-dominant part of the primordial density perturbations. This fact was used in Ref. 15 1 
to infer upper bounds on the duration of inflation. These bounds are compatible with ad- 
equate inflation to solve the horizon and flatness problems, but raise the possibility that 
the total duration of inflation might be observable. This possibility goes against a com- 



monly held view that inflation erases the memory of anything which occurred previously, 
and hence increasing its duration beyond the minimum needed to solve the horizon and flat- 
ness problems can produce no observable effect. However, contrary indications to this view 
had previously been published in the form of arguments that inflation cannot be eternal to 
the past[l6|, llTJ, although these arguments were based on general considerations which do 



18| has recently suggested a 



not make speciflc predictions of observable effects. Winitzki 

model in which inffaton fleld fluctuations can produce violations of the null energy condition 

and possible effects of the total inflationary expansion. 



The purpose of the present paper is to re-examine and improve the analysis in Ref. |l5 |. 
In Sect, ini we discuss a kinematic model which makes no explicit reference to the inflaton 
field, but examines the gravitational effects of the stress tensor fluctuation upon timelike 
geodesies. In Sect. Illlt we give a detailed treatment of a dynamical model in which the 
stress tensor fluctuations alter the dynamics of a scalar inflaton fleld. In both models, a 
correction to the power spectrum of density fluctuations is computed. Section |IV] discusses 
the implications of our results and some associated conceptual issues, especially the role of 
transplanckian modes. Our analysis is summarized in Sect. El 

Before concluding this section we should mention some conventions. A hat is used to 
denote the spatial Fourier transform of any fleld A(t,x.) 

1 



A{t,k)^—--Jd'xe'''''A{t,^). (1) 

We represent the power spectrum of A by the symbol VA{k,t), which is deflned as follows 
from the correlator of two A flelds: 

(i(t, k)i(t, k')) = VA{k, t) X ^^±^ . (2) 

In (the usual) cases for which A{t, k) is time independent we drop time from the argument 
list of the power spectrum, as in P^(fc). We consider the loop counting parameter of quan- 
tum gravity to be k^ = 16ttG. Our curvature tensors follow the Landau-Lifshitz spacelike 
convention, which is also the Misner-Thorne- Wheeler (+++) convention. 
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A very important point for understanding our analysis and results is that we normalize the 
FRW scale factor to unity at the end of inflation, rather than at the current time. Also 
note that we use the subscript "0" sometimes to signify "background" and sometimes to 
denote that the subscripted quantity is evaluated at the beginning of inflation. So to is the 
time at which inflation begins, rather than the current time as in much of the literature 
on cosmology. We indicate the current time by the subscript "now", so the wave number 
k = 2it/X is measured in units of the comoving distance at the end of inflation, and it can 
be expressed in terms of the current wave number /cnow = 27r/Anow through the relation 

rC '3'now'^now 

II. THE KINEMATIC MODEL REVISITED 



Here we will review and modify a model flrst presented in Ref. 15|. The point is to give 
a simple computation of the extra part of the power spectrum of energy density fluctuations 
due to a conformally invariant quantum fleld. (See Fig. [1] for the relation between our 
contribution and the usual tree order result.) A rigorous derivation involves solving the 
coupled, linearized inflaton-graviton equations with the conformal stress tensor as a source. 
We will do that in Sect. IIIII Here we avoid any mention of the inflaton fleld and we require 
only the background metric 

ds^ = -dt^ + a^{t) 5ijdx'dx^ = a^{r]) {-drf + 5ijdx'dx^) , (4) 

where t is the comoving time, 1] is the conformal time, and a is the scale factor. 

What we do instead is to assume that the stress energy consists of a perfect fluid with 
energy density p(t, x), pressure p(t,x) = wp(t,x) (with constant equation of state w) and 
4- velocity 'U^(t,x), in co-moving coordinates such that u'^d^ = d/dt. Then we use energy 
conservation, 

P + {p + p)e = o, (5) 
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FIG. 1: Various contributions to the power spectrum of primordial perturbations. Wavy lines 
stand for graviton-inflaton fields and solid lines denote conformal fields. The leftmost diagram 
represents the tree order contribution which is usually reported. The center diagram gives the 
one loop contribution from conformal matter which is the subject of our analysis. The rightmost 
diagram represents the (unobserved) term which is neglected by subtracting off the expectation 
value of the conformal stress tensor. 

with p = dp/dt, to infer the perturbed energy density Sp{t, x) by perturbing the expansion 

The key to simplifying the computation is deriving the perturbed expansion 6d{t,x.) from 
the Raychaudhuri equation 

u>'d,9 = -R^^u^u-' - ^9^ - a.^a^'' + u'^^'u^^ + {u^u^),^ . (7) 

We shall drop the shear a^'^ , the vorticity uj^'^ , and the acceleration (■u'^.^-u'^);^, at which 
point one can obtain 50{t^ x) from SRfj_u{t, x) = |K^T™°^(t, x), for the particular part of the 
total perturbation that concerns us. (See Fig. [H) Here T^°"^ denotes the stress tensor of the 
conformal field. 

This makes for a wonderfully simple analysis in which we need never consider the per- 
turbed inflaton field or components of the perturbed metric. Unfortunately, it isn't correct, 
as we will see in Sect. Illli Ignoring a'^'^ and w^'" is valid at linearized order for single-scalar 
inflation, but the acceleration term contributes at linearized order and that spoils the simple 
relation between 66{t,x.) and the conformal stress tensor. So the result we shall derive in 
this section is off by an important factor of {k/H)^, but it does depend correctly on the 
initial time. 



Conformal invariance allows the stress tensor correlation function 

to be written in terms of the flat space stress tensor correlation function. 

CS(^, ^') = a~\v) a-'W) CilU^, x') . (9) 

Here the components of C^J^^(x, x') are understood to be in the second set of coordinates 
in Eq. (|1]). Although the conformal anomaly term in (r^y{x)) breaks conformal symmetry, 
this term cancels out of the correlation function, Eq. ([H]). 

The fiat spacetime energy density correlation function of the conformal field is S{Ari,r), 
where Arj = rj — rj' and r = |x — x'|. The expansion correlation function can be expressed 
in terms of S{Ari, r) as 

(5e(r^i,x) 5eiv2,^')) = -n' a'^r^,) a'^r^,) T ^ T -% £{Av,r) . (10) 

4 Jrto a{r]) Jrio a{r]') 

For the case of the electromagnetic field, 

and the expression for the conformal scalar field case is identical except for an additional 



factor of 1/12. (Note that this expression corrects an error in Eq. (39) of Ref. 15|.) This 



ll|. 



expression is ultraviolet finite in spite of being one loop 

After reheating, the expansion fluctuations cause differential redshifting and consequent 
density fluctuations, in accordance with the conservation law Eq. ([S]) for a perfect fluic . 



The fluid flow approach to density perturbations has been discussed by several authors 19- 



22|. Let 6p = 6p/p be the fractional density fluctuation at conformal time rj = rjs, the 
last scattering surface. Its spatial correlation function is given in terms of the expansion 
correlation function by 

{6,{vs, x) 6,{vs, x')) = (1 + wf r -^ r -^ {6e{r^u x) 6e{^,, x')) . (12) 

Here reheating occurs at i] = rjr and the equation of state after reheating is p = w p. 



Let 

Fo{r) = {6e{r]r,^)6e{r]r,^')) (13) 

be the variance of the expansion at the end of inflation. In many cases, the dominant 
contribution to the density fluctuations arises from effects occurring during inflation. Then 
contributions to the expansion correlation function coming from stress tensor fluctuations 
after reheating may be neglected. However, we still need to account for the evolution of 
66 after reheating. If we ignore the effects of classical density perturbations and pressure 



gradients, as well as the quantum stress tensor, then 69 satisfies [See Eq. (13) in Ref. |l5|.] 

where ^o = 3a/a is the unperturbed Robertson- Walker expansion. The solution of this 
equation can be written as 

S^iv) = ^ , V>Vr, (15) 

where we have set the scale factor at reheating to unity, a(?7r) = 1- In the approximation 
where we consider only stress tensor fluctuations during inflation, after reheating we have 

{6e{ri,,x)6e{ri2,x')) = a'' (ri,) a-"" (m) Fo(r) . (16) 

The density perturbation correlation function now becomes 



((5„(r/„x)5„(r/„x'))^(l + w;)'Fo(r) 
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'?" dT]i 



(17) 



The time integral now depends only upon the form of the scale factor between reheating 
and last scattering. We consider a model in which the inflation is described by a de Sitter 
metric, 

a{v) = -j^, V<Vr, (18) 

and the subsequent period is radiation dominated {w = 1/3), 

a{ri) = Ht] + 2. (19) 



Here H is the Hubble parameter of de Sitter space, and both a{ri) and da/drj are continuous 

a.t 1] = Tjr = —l/H. Then 

fVs dm p^ dm 1 , , 

' ' A ~. (20) 



The last step follows because a{ris) ^ 1. Now we obtain 

((5,(r/„x)(5,(r/„x')) = ^Fo(r). (21) 

We will next take spatial Fourier transforms and write, for example, 

£U^V,k) = -^ld'xe'''-^^£U^r],r). (22) 



The explicit form of Eemi^f]-, k) is 



- k"^ sin{k Ai]) k^ [^ 

^^-^^^' ^) = ' 960.^ Ar, = -960^^ Jo ^^ '°'^^" ^^^ " ^'^^ 

The Fourier transform of Fo(r) is 

AW = i.W"4Lr^5(A,,,). (24) 

4 J,,o a(r7) 7% a(77') 

We next evaluate the integrals in this expression, using the second form in Eq. 023]), to 
find in the limit that k\T]o\ ^ 1, 

^»w-^0"»i' + i 1*1^ + ■■■)■ (^«' 

Note that the magnitude of -Fo(^) grows as Ir^ol — t- oo for fixed k. The growth was found in 
Ref. (l5|, but there only the jr/op term appeared. The reason for this is that in Ref. jlSJ, the 
calculations were done in coordinate space until the last step, where a term in -Fo(r) oc 1/r^ 
was Fourier transformed into a term in Fo{k) oc k^. However, this procedure is not sensitive 
to the possibility of delta-function terms in Fo{r). The leading term in Eq. (!25|) arises from 
just such a term, one proportional to V^5(x — x'). 

Let Psp{k) denote the spatial Fourier transform of (5p(?7s,x) (5p(?7s,x')), which is related 
to the power spectrum by Vs {k) = Ank^ Ps {k). From Eqs. fl2T|) and fl25|) . we find 

Ps^ik) ^ ^/ci^.e-^— ')(5,(,.,x)5,(,.,xO) - ^^ [-\krjonl\k,on...) . 

(26) 



Note that if we take the k'^ |?7o|^ term seriously for all k, then its spatial Fourier transform 
contributes a term proportional to V^(5(x — x') in {6p{ris,^) 5p(?7s,x')). To the extent that 
measurements are made in position space by comparing proxies for 6p{t, x) at x 7^ x', this 
term will not contribute. Here we assume that we may ignore this term and retain only the 
nonlocal t]q effect. Recalling the definition Eq. ([2]), our result for the one loop contribution 
to the 6p power spectrum from conformal matter is then 



Apart from numerical factors, this is equivalent to Eq. (48) in Ref. 15|. We will see in the 
next section that the prefactor of k'^ should really be if ^, however, that still leaves the result 
very strongly biased toward short wavelength perturbations and far from scale invariant. 
The possible implications will be discussed in Sect. [IVl 

III. A DYNAMICAL MODEL 

The computation of [Vsp{k)]coni we have just completed involved three basic steps: 
• Inferring the post-inflationary density contrast from the perturbed expansion 

6p{t,^) = -{l + w) f dt'Se{t,x.). (28) 



• Approximating the post- inflationary Raychaudhuri equation as flT4l) . so that the den- 
sity contrast during radiation domination becomes 

^p(^,x)^-^^5^(t„x). (29) 

• Computing the perturbed expansion which is accumulated during inflation by using the 
Raychaudhuri equation ([71), under the assumption that the acceleration term {u^.^^u'^) -^ 
makes no contribution at linearized order. 

Equation fl28|) is exact. While Eq. ( 129!) is certainly not exact, it does represent a reasonable 
approximation when Fourier transformed and restricted to super-horizon modes. The prob- 
lematic step is ignoring the acceleration term to compute the 66 [tr, x) induced by conformal 



matter fluctuations during inflation. It turns out that the acceleration term depends linearly 
upon the inflaton perturbation, and we must study the coupled gravity-inflaton system to 
get a reliable result for 69{tr,x.). Having done this, we use Eqs. f l28|) and fl29l) as before to 
compute [Vsp{k)]coni- 

In this section we study the coupling of a single scalar inflaton field with the gravita- 
tional perturbations of a spatially flat Robertson- Walker spacetime. These perturbations 
are in turn driven by the fluctuations of the stress tensor of a conformal quantum field. The 
unperturbed metric is of the form in Eq. (jl]). During inflation, this metric will be approx- 
imately that of de Sitter spacetime, although with a slowly varying Hubble parameter H. 
We assume that the unperturbed metric satisfies Einstein's equations with the stress tensor 
of a spatially homogeneous inflaton field <fo{t) as the source. Let the inflaton be self-coupled 
by a potential V{!f). Then the Einstein equations become 

5H' = ^k' (^v^o + K)) , (30) 



and 



-2H- 3H' = \^' {^^^l -Vo). (31) 



Here dots again denote derivatives with respect to t, H = a/ a and Vq = V{(po). The scalar 
field equation is 

^ -.dj^^g'^'-d,^) + V'iv) = , (32) 



-9 



which becomes 



y3o + ^HSipo + V^o' = . (33) 

To the extent that Vq is not constant, the unperturbed spacetime will not be exactly de 
Sitter space. 

A. Coupled Equations for Inflaton and Metric Perturbations 

We next wish to consider linear perturbations of this spacetime in a gauge in which 

9u = -l, (34) 

10 



so the perturbed metric may be written as 



ds = -dt + 2a{t)hu{t, y.)dtdx'- + a (t) 6ij + hij{t, x 



dx^dx-' . 



(35) 



It is convenient to define the conformally transformed, spatial metric, 



9ij = ^ij + h 



ij ■ 



(36) 



The determinant of the full metric can be broken up into three factors. 



g = a'^ X det{g) x 1 + huhtjg'^ = a^ l + h + 0{h^) 



(37) 



where h = 5^^ hij is the trace of the metric perturbation. 

In addition to the metric perturbation, the infiaton field will have inhomogeneous per- 
turbations: 

ip{t,^) = ipo{t) + 6ip{t,^). (38) 

The scalar field perturbations will satisfy 



5(^ + 3H5ip + -ipoh + VQ5ip = , 

which follows from the expansion of Eq. (132|) to first order both in S^p and in h. 
The Einstein equations may be written as 



R^, = ^-^ fe^^ - ^-g.^g^'X''^^' 



(39) 



(40) 



where the total stress tensor is the sum of contributions from the infiaton field and the 
conformal quantum field: 



T-itotal T-iinfl I T-iconf 



(41) 



We focus here on the time-time component of the Einstein equation. The first order expan- 
sion of Rtt is 



Ru = -3H -3H^-^(h + 2Hh) + -[hu,^ + Hhu,i) + ©(/i^) . 



(42) 



The infiaton stress tensor satisfies 

1 



^f - o^.^^'"'rr = d,^d.^ + g,.V{v) . 



(43) 
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The first order expansion of tlie tt component of this expression is 



dtifdt'f - V{ip) = 00 - V{(fo) + 2(^o5ip - V'{ipo)Sip 



Thus the equation for the first order metric perturbation can be written as 



(44) 



1 1 K^ 

- -(/i + 2Hh) + -{hti^i + HhtL^ = y [20o<^0 - V'{ipo)S^ + U 



where we define 



U = T, 



conf 

tt 5 



(45) 



(46) 



the energy density of the conformal field in the comoving frame. 
Define the normal vector to the surfaces of constant ip by 

The first order expansion of the spatial components of this vector is 

i hti di5ip 



(47) 



We can impose the gauge condition 



«'^ = '5f, 



(48) 



(49) 



from which the condition gu = —1 follows. In addition, this leads to a relation between 6ip 
and hti to first order: 

hti{t,yi) = -di 



5ip{t,x) 



_a{t)ipo{t) 
This relation allows us to eliminate the hu terms in Eq. (H5|) and write 

1 d 



(50) 



1 

— / 
2 



a? dt KipQ 



^^o5^ + ^K^V^S^ - ^^ ( ^ ) - 7^^ {a'^h) = -f^ 



U 



(51) 



2a'^dt ^" '' i 

Equations (l39ll and ( ISTl) form a pair of coupled second order equations for the metric and 
scalar field perturbations. These equations may be rewritten by expressing y^o, Vq, and their 
derivatives in terms of the Hubble parameter Hit) and its derivatives. The sum of Eqs. ( !30|) 
and (13T1) leads to 

K(pQ = 2V^ , (52) 
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and their difference leads to 



From these relations, we find 



/tVo = 2H + QH'^ 



(53) 



H 



and 



kv;, 



Kipo 



H 



-H 



6H\^-H 



Kifo 



V" 



H 



ij2 



-H 2(-if)f 



(54) 



H H^ 3HH ■ , , 

- + —J7-. - —^ - ^H . (55) 



_jj 2H 4if2 2H 

Note that the homogeneous form (setting [/ = 0) of Eqs. ( 139|) and ( 15T|) has a solution 
when 6(f = 0o and h = 6H. We can reduce the order of the system by scaling out this 
solution and defining 



B(x) 



6ip{x) 



(56) 



We can now rewrite Eq. ( l39l) as 



.V^o 



h = 6HB -2{2— + 3H]B-2B. 
V 00 J 

This allows us to eliminate h from Eq. (ISTI) . and write 

O B = -k'^U, 
2 

where O is the operator defined by 
O 



(57) 



(58) 



d^+ Uh + 2^] d, + 4H + 6H' + 4H^ -2^ + 2^ -\ 
\ 00/ 00 00 00 a , 



(59) 



Equation (158|) is a third-order equation which we will solve using a retarded Green's function. 
It will be convenient to take a spatial Fourier transform, and define the operator 

O, = d^ + (5H + 2^] d, + 4H + 6H' + 4H^ -2^ + 2^ + ^ . (60) 

V V^o/ V^o V^o V^o a^ 

Let G(t, t', /c) be the retarded Green's function of this operator, which satisfies the equation 



OkG{t,t\k) = 5{t-t') 



(61) 
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with the boundary condition 



G{t,t',k) = if t<t'. 



(62) 



Let \l/i and \l/2 be two hnearly independent solutions of the homogeneous equation 



0fc*(t,fc) = O. 



(63) 



The Green's function may be expressed as 



1 



G{t, t', k) = — ^ [vl/i(t<, fc)vl/2(t>, k) - ^^{t, k)^2{t', k)] , 
where t< and t> are the lesser and the greater, respectively, of t and f, and 

W{t, k) = ^i(t, k)^2{t, k) - ^i(t, fc)^2(t, fc) 



(64) 



(65) 



is the Wronskian. 

The homogeneous solutions \&i are difficult to obtain in general. However, if we make a 
"slow roll" approximation in which time derivatives of if, and hence of ^Jq and of Vq, are 
assumed to be small, then we have approximately 



Ok ^ d"^ + 5H dt + QH^ 



In this approximation, the solutions of Eq. f l63|) are 



i2 ■ 



(66) 



vl>^(t,A;) = a-2(t)e^'^4'^*^""(*^) 



(67) 



and ^1^2(^5 k) = "^{{t, k). Here to is an arbitrary constant. Now the Wronskian is 

2ik 



W{t,k) 



a^{t) 



(68) 



and the retarded Green's function may be written as 



G{t,t',k) = -^^^ sm 



A;a2(t) 



* dti 



(69) 



for t > t'. 
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B. Inflaton Field Fluctuations 

In this subsection, we wish to calculate the fluctuations in Lp which are driven by the 
stress tensor fluctuations of the conformal field. Let G{x^ x') be a coordinate space Green's 
function for the operator (9, which satisfies 

O G(x, x') = 6{x - x') = 6{t - t') (5(x - x') . (70) 

A particular solution of Eq. S5E\i can be written as 

1 r 

B{t, x) = -K^ / d^xi G{x, xi) U{xi) . (71) 

We now treat U and hence B as fluctuating fields, and write the correlation function for 
B as 

{B{t, x)5(t, x')) = jK^ f d^xid^X2 G{x, xi)G{x', X2) {U{xi)U{x2)) . (72) 

Next we convert from comoving to conformal time, using drj = dt/a{t), and use the rela- 
tion between the (comoving) energy density in Robertson- Walker spacetime to that in flat 
spacetime, 

{U{xi)U{x2)) = a-'(r/i) a-\r]2) S{Ar], r) , (73) 

where S{Ari,r) is the flat spacetime energy density correlation function, with r = |xi — X2I. 
This leads to 



l,.4^^^^^^^/^ p dr]i p' dr]2 
X [G(77, r/i, X - xO G{i, r^^, x' - x^) ^(Ar^, r)] . (74) 



(9,5(r/,x)9,i?(V,x')) = -K'a{T^)a{i) f 4N f' ^t^J d'xj'x2 



Here the boundary condition drfB{ri,x.) = at rj = rjQ has been imposed. 
Next we take spatial Fourier transforms, and define 

{d,B{7], x) d,,Bir]', x')) = / d'k e^'^-^^'^'^ (9,5 d,,B)k , (75) 

and analogous transforms of G(77, 771, x — xi) and of £{Ar], r). Then we may write 

4 J no a-^{r]i) Jrio a-^{r]2) 

G{7],7],,k)G{r]',r]2,k)S{Ar],k)] . (76) 

15 



The correlation function for drjB may be integrated to yield the mean squared inflaton 
fluctuation at the end of inflation, tj = rj^. \i B = Q aX rj = tjq, then 

{B\r]r))k = r dr] T dr]'{dr,B d^,B)k . (77) 

Jrio Jr]o 

The form for the Green's function from the slow roll approximation, Eq. fl69|) . may be 
expressed as 

G'(r;,V,A:) = -|^ sm[kiv - v')] ■ (78) 

With this form, and the de Sitter space scale factor, Eq. f lTS]) we find 

{B'^{r]r))k = —rrr / drjri / dr]' rj' / drji sm.[k{ri - rji)] 

X r dr]2 sm[k{r]' - r^s)] S{Ar], k) . (79) 

Jvo 

It is interesting to compare this with the result Eq. ( ITOl) of the kinematical model. After first 
horizon crossing, 69(1], k) ~ H^B{rj,'k) [See Eq. (l86l) below.], so the big difference between 
Eqs. (179|l and ( ITOj) is the extra factors of 

—— I drji sm[k{ri-rii)] x — - / dr]2 sin [A; (77' -772)] . (80) 

fc Jr}o rC Jr}o 

These terms describe how stress tensor fluctuations from very early times are communicated 
by the inflaton field to the late time geometry, and they effectively introduce a factor of 
{H/kY X {H/ky = {H/kY to the result of the kinematical model. 

Finally, we may use the form of £{Ar],k) given in Eq. (!23|l . and perform the integra- 
tions in Eq. (179|) using the algebraic computer program Mathematica. The result is rather 
complicated, but in the limit of large |?7o|, it becomes 

^4 0-2 |„ 13 ^4 tt2 2 

\ Kirj/k i228807r4 ISSGOOtt'^A; ^ ^ 

As with our result Eq. fl25l) for Fo{k), the leading contribution for large |?7o| corresponds to 
a term which is ultralocal in position space, in this case proportional to (5(x — x'). 

In the subsequent analysis, we will also encounter expectation values of quadratic forms 
involving time derivatives of B, such as (BB) and (-B^). However, one may check that all 
of these terms are at most proportional to |?7o| and hence sub-dominant compared to (5^) 
in the limit of large \rio\. 

16 



1. Density fluctuations from the conservation law 

As discussed at the beginning of this section, we can compute the density contrast using 
Eqs. (128!) and ( 129|) once we have 69(tr,yi), the perturbed expansion at the end of inflation. 
The expansion 6 can be obtained from its definition 

e ^ u'^.^^ = -^d,{V^ u^) . (82) 

In our gauge, Eq. ( 149|) . this becomes 

0(t,x) = |-ln(v/=^). (83) 



dt 
We may use Eq. (!37|) to write 

9 = 3H+l^ \n{g) + 1|, ln[l + huht.r] = 3H + h + 0{h^) . (84) 



2dt '^' 2dt 



2 



Recall that Oq = 3H is the expansion of the comoving geodesies in Robert son- Walker space- 
time, so the first order perturbation of the expansion is 

50(t,x) = i/i(t,x). (85) 

One infers h{t,x) from B{t,x.) using Eq. (157|) . Because we only need it at the end of 
inflation the B and B terms can be dropped, and we can use the radiation domination 
result H = —2H^ to conclude 

6e{tr, x) ^ -QH'^B{tr, x) . (86) 

Hence Eqs. (1251) and (12^ give the following expression for the density contrast during radi- 
ation domination: 

5,(t,x)^-4i75(t,,x). (87) 

This should be valid when Fourier transformed and restricted to super-horizon modes. 

To find the power spectrum we first compute the spatial Fourier transform of the 6p 
correlator using Eqs. (!87|) and (ISTj) 

^'"^^^ = (^/rf'^e^'(^--') (5,(77., x)5,(r^„x')) ^ -j^^:^ [-\kr^,\^ + -\kr^o? + ■ ■ ) ■ 

(88) 
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Multiplying by Airk^ gives the power spectrum. As for the kinematic model (126|) . we assume 
that we may drop the |?7oP term which is ultralocal and presumably not part of the observed 
power spectrum. This leaves us with 



VsAk) 



(89) 



conf 2400 7r2 \Ha{to)J 
Thus the dynamical model also produces a non-scale invariant spectrum biased toward the 
blue end of the spectrum, although less so than in the case of the kinematic model. 



2. Density fluctuations from the Sachs- Wolfe effect 



An alternative approach to calculate density or temperature fluctuations is to study the 
effects of metric perturbations on the redshifts of photons, as was first done by Sachs and 



Wolfe 



24| . Equation (39) of their paper may be expressed as 



AT 



dt 



e'hu^t{x) - -e'e^ hij^t{x) 



(90) 



This formula gives the differential redshift, and hence temperature fluctuation, of a photon 
propagating from t = tr to t = tg along a null geodesic in the direction of the unit vector 
e\ The integrand is understood to be evaluated along the unperturbed null geodesic. In 
contrast to the previous discussion, we now need expressions for the individual components 
of the spatial metric perturbation, h^. For this purpose, it is convenient to express the 
scalar part of the perturbed metric, Eq. fl35|) as 



ds' 



dt^ - 2B^idx'dt + 0^(1 - 2ip)Sij - 2E, 



V 



KAjtXj KajJLj • 



(91) 



Here we follow the notation of Mukhanov 



6|, as modified in Ref. |25|. As before, hu is given 



by Eq. ( l50l) . The spatial components of scalar perturbations are 



hij{t,x.) = -26ijilj{t,:>i) 



■did,E{t,^). 



(92) 



a2(t) ^''^' 

The quantities which appear in the integrand of Eq. (!90i) may be written in terms of B, ip 
and E as 



e'hti^tix) 



e-V 



[B - HB) 



(93) 
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and 

- h'& h^^^tix) = ^ + (^) {E - 2HE) . (94) 

Expressions for ip and for E — 2HE, Eqs. (1A23|) and (]A25p . respectively, are derived in the 



Appendix. First we note that E — 2HE contains two types of terms, those which involve B 
and B, and those which depend upon U or U evaluated at the same time as E — 2HE. Both 
of these types of terms will give a sub-dominant contribution, which is either independent 
of |?7o| or small compared to the \7]o\^ and |?7oP terms. The same comment applies to all 
terms in ip, except for the H B term. Note that in coordinate space, Eqs. ( 1A23I) and ( 1A25I1 



contain the non-local operator 1/V^, which is non-local in space only, not in time. In any 
case, calculations are best done in Fourier space, where 1/V^ is replaced by —l/k"^. 
Thus we may take 

ijj^HB. (95) 

If we drop the B term in Eq. ( l93l) . and use the fact that 

d d /e-V\ 

Jt-m-[—) ("«' 

is the total derivative along our null geodesic, we may write 

In the last step, we used the fact that the dominant contribution will come from the lower 
limit of the integral. If we recall that here 

we again obtain Eq. (IHTj) . 

Density perturbations are often treated using the gauge invariant potentials, which is 
yet another possible approach. However, both the fluid flow approach and the Sachs- Wolfe 
formula, Eq. (MI]) , are themselves gauge invariant and somewhat simpler for our purposes 
than the gauge invariant potentials. 
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IV. DENSITY PERTURBATIONS FROM QUANTUM STRESS TENSOR FLUC- 
TUATIONS 

A. Possible Constraints on the Duration of Inflation 

Let us now discuss the possible physical implications of the conformal matter contribution 

to the power spectrum. Recalling that k^ = IGtcG, and that current wave numbers k^^ow 

correspond to k = {ak)now, our result fISI?]) can be expressed as 

r . 8G^({ak^y 

^''"Jconf-^ 75 [ aoH J ■ ^^^^ 

This is not scale invariant, and it is associated with highly non-Gaussian fluctuations. In 
contrast, observations of large scale structure and the cosmic microwave background radi- 
ation are consistent with the primordial perturbation spectrum being approximately scale 
invariant and Gaussian 71] 

/ , \ -0.037±0.012 

P,(t„„,) « (2.4411-1) X 10- (;^35f=-r) . (100) 

The primordial curvature and density contrast power spectra are related by V-n = j^Vsp 
6|.) Note that the weak scale dependence of the observed power spectrum f llOOp is actually 
in the opposite (red) sense to the massive blue tilt we predict from conformal matter. Hence 
the contribution from conformal matter can only represent a tiny part of the total power 
spectrum. Because our result Eq. f l99|) grows like l/a^(to) as the start of inflation is pushed 
back to earlier and earlier times, one can derive a bound on the duration of inflation by 
requiring that Eq. (1^ is small enough to not affect the measured result, Eq. fllOOp . 

It will facilitate the discussion to recall some reasonably generic predictions of single- 
scalar inflation in the slow roll approximation. The tree order results for the scalar and 



tensor power spectra are 22] 



Vnik^o^)]^ ^ ^^'^y , [P.(A:now)l^ ^-GH^t^), (101) 

J tree 7re(tfc) L -I tree -ji 

where e(t) = —H/H^, and tk is the time of flrst horizon crossing, 

flnow/^now = a{tk)H{tk) ■ (102) 
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The absence of much scale dependence in the observed result, Eq. fllOOl) . is explained by H{t) 
being approximately constant during inflation. (This is why our de Sitter approximation 
of Sect, mil was well motivated.) Of course a nearly constant H{t) makes the slow roll 
parameter e(t) = —H/H^ close to zero. The enhancement of the scalar power spectrum by 
l/e(tfc) explains why it has been observed, while the tensor contribution has so far not been 



resolved. At 95% confidence the bound on their ratio is 



i 



Pft(0.002 Mpc-i) 
^(a002Mpc^ 



^ ^ r.rr. 7.-1 < 0-22 ■ (103) 



With the theoretical results Eq. (1101 p and the scalar observation Eq. (llOOp . this implies a 
bound on the inflationary Hubble parameter 

GH^ = — X r X ^7^(0.002 Mpc"^) < 10"^° . (104) 

16 

Note that our one loop contribution Eq. ( l99l) is suppressed by GH^e relative to the tree 
effect Eq. (llOip . It can only become observable when inflation begins at such an early time 
that these factors are canceled by the square of the physical wave number in Hubble units. 

The bound we get on to derives from requiring the predicted contribution from conformal 
matter Eq. ( l99l) to be smaller than the observed result Eq. (llOOp for the smallest wave 
number k^ow for which data exists. We take k^ow ~ lO^^^cm^^ ^ 2 x 10"^^ GeV, which 
corresponds to structures of about 2 Mpc in physical size, or about 5 arcminutes of angular 



scale 
that 



23|. Let Tr stand for the reheat temperature, and let us assume efficient reheating so 



i/^ ^ SvrGT^ ^ Svr X 10^° GeV^ (^^^^) . (105) 

The Universe has expanded by about a factor of 10^ since the last scattering time tg (when 
the temperature was Ts ^ 1 eV) and recall that we normalize the scale factor to unity at 
the end of inflation, hence 

a_ - 10^ aiQ - 10^ g - 10- (y^t^) • (106) 
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For these parameters our result Eq. (p9|) implies 

Requiring this conformal contribution to be less (by a factor of ten, say) than the observed 
spectrum Eq. fllOOp gives 

i<io«(i2;^)'. ,108, 

CtO \ Tr J 

Recall that sufficient inflation to solve the horizon and flatness problems requires 1/ao ^ 10^^, 
so Eq. f llOSp allows more than enough inflation for this purpose. Note that this bound is, 



apart from being improved by a factor of 10'^, equivalent to Eq. (92) in Ref. 15|. However, 
the latter result was derived using an overly simplifled dynamical model which did not fully 
account for the coupling between the inflaton fleld and the perturbations of the spacetime 
geometry. 

B. The Transplanckian Issue 

We now turn to some of the conceptual issues which are raised by the calculations de- 
scribed in the previous sections. One of these concerns the use of transplanckian modes, that 



is, modes whose physical wavelengths are less than the Planck length of ip = yG ~ 10^'^^ cm, 
as measured by an observer at the start of inflation. The estimates in the previous subsection 
dealt with perturbations with a present wavelength on the order of Anow ~ 10^^ cm k. 10^^£p. 
The physical wavelength of this mode at the beginning of inflation is, 

Ao = ctoA = ao X -^^^ . (109) 

'^now 

If we assume that oq is at the bound (llOSp — which means conformal matter contributes 
10% of the measured power spectrum at the smallest observed scales — then the initial 
wavelength is 

For all reasonable values of Tr, this is far below the Planck length. 
This raises two questions: 

22 



• Is it valid to extrapolate low energy dynamics such as electromagnetism to trans- 
planckian scales? 

• Is it valid to apply perturbation theory for transplanckian modes? 

No one knows what dynamical principles might apply at Planck scales, but it is of course 
acceptable to carry out a study, as we have done, based on the explicitly stated assumption 
that they are unchanged. What doesn't seem alright is employing perturbation theory. One 
must not be mislead by the fact that the tree order effect ~ k^H"^ /e is small; the series is 
an expansion in powers of the large parameter [Kk/aoY, 

pu*)~-'W^+«.(^)%«.(^)%...|. (Ill) 

If one makes the usual assumption that the pure numbers a^ are of order one then the only 
way of making the one loop term comparable to the tree order result must also make the 
two loop and higher terms comparable. The conclusion seems unavoidable that perturbation 
theory must break down, for mode fc, as the initial time is pushed back to the point for which 
K,k/aQ ~ 1. We do not possess a nonperturbative computational technique so what actually 
lappens at earlier times is a matter of conjecture and lively debate within the community 



26 



28|. 



One view is based on the observation that the far ultraviolet contains so many modes 
that even very small deviations from quiescence in each of them must produce enormous 
fluctuations that would invalidate semi-classical inflation. Hence it must be, the argument 
goes, that a nonperturbative resummation of loop corrections such as Eq. f illip exhibits no 
large effect, even for very early initial times. For each wave number k there would be a time 
Tfc such that nk/a{Tk) ^ 1, after which our perturbative treatment is valid. As long as to 
comes after T^, making to smaller causes the one loop effect to grow as we predict, with 
higher loop contributions still negligble. But if to is pushed before T^ then the higher loop 
corrections become important and the whole series approaches a constant. If this view is 
correct then, for to < T^, one could only employ the perturbative treatment of this paper by 
starting the evolution of i?(t, k) at t = T^., not at t = to- And the correct initial condition 
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would be something close to quiescence at t = T^. This is a nonlocal initial condition, but 
then quantum effects typically are nonlocal. 

A different view is motivated by the similarity of these issues to those which arise in black 



hole physics. The original derivation of the Hawking effect 29| assumes free quantum field 
theory on a fixed background spacetime and requires transplanckian modes. This derivation 
is analogous to our treatment in the previous sections. It is possible to reproduce the 
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31| . but only by postulating a 



Hawking effect without the use of transplanckian modes 
non-linear dispersion relation, which breaks local Lorentz symmetry. 

If there is a new physical principle which avoids transplackian modes, then economy of 
thought would suggest that it should be the same principle for both black hole physics 
and for cosmology. Ideally, one might hope for an experimental or observational test of 
transplanckian physics. The power spectrum which we have derived using transplanckian 
modes has the potential to provide such a test. If inflation lasted just slightly less than the 
amount given by Eq. f llOSp . then the model described above predicts a non-scale invariant and 
non-Gaussian component in the cosmic microwave background which might be detectable. 

C. Relation to Weinberg's Theorem 

Neither the comoving wave number k nor scale factor a{t) are physical, only their ratio, 
k/a{t). Of course ratios of the scale factor at different times are also physical. Because 
we normalize the scale factor to one at the end of inflation, the various factors of k in our 
results must really be interpreted as the physical wave number at the end of inflation, k/ar- 
Therefore, the kinematic model estimate of ['P<5p(/i;)]conf ~ K^k^/[H'^al] seemes to suggest 
a one loop correction to the power spectrum which not only violates scale invariance by 
the factor (k/ar)^ but also grows at late times like (0^/00)^. Such growth would contradict 



a bound of at most logarithmic growth established by Weinberg 3^. (Weinberg's result 
was derived for minimally coupled scalars but it can easily be extended to conformally 



coupled particles 33|].) In fact one can see from the dynamical model of Sect. IIIII that 



24 



there is no growth at late times; what happens instead is that the principal effect arises 
from fluctuations near the time to when the interaction is turned on, after which it rapidly 
approaches a constant. By comparing our one loop correction with the usual tree order 
result 



k'H 



Vsik)_^^ versus Vsik)_^K'H-x-^, (112) 



K'k' 



conf 



tree £ 

It will be seen that our contribution consists of the tree result (without the inverse of 
e = —H/H'^), multiplied by a typical one loop correction of the square of k, times the 
mode's physical energy at the initial time. Later times contribute far less because the 
mode's physical energy redshifts so rapidly. There is no mystery about why the effect can 
be large at very early times because the mode is transplanckian then and should induce 
large gravitational effects. Of course this again raises concerns about using perturbation 
theory and low energy dynamics. What Weinberg considered was quantum loop effects from 
the "safe" regime of late times during which perturbative general relativity must be valid. 
Our results are in perfect agreement with his bound; indeed, they fail to show even the 
logarithmic growth allowed for by the bound and achieved by nonconformal matter. 

V. SUMMARY AND DISCUSSION 

In this paper we have evaluated the extra contribution to inflationary density perturba- 
tions from the quantum stress tensor fluctuations of a conformal field such as the photon. 
This was done in a simple, kinematical model and then in a more accurate, but much more 
complicated, dynamical model. Our main result is that the power spectrum of the energy 
density at wave number k goes like the tree order result (without enhancement by 1/e) 
times {E{to)/Mp\)'^, where E{to) = k/aito) is the mode's physical energy at the beginning 
of inflation and Mpi is the Planck mass. If a perturbative computation such as this could 
be trusted to arbitrarily early times, the absence of such a massive blue tilt in the observed 
power spectrum would seem to imply a bound on the total duration of inflation. This 
constraint allows enough inflation to solve the horizon and flatness problems. 
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Our result derives from ver y ea rly times and rapidly approaches a constant, so it does not 



contradict Weinberg's bound [32 



33| , on the maximum possible growth of quantum correc- 
tions at late times. However, our result does involve a problematic extrapolation of known 
physical laws to the transplanckian regime, and the even more problematic assumption that 
perturbation theory can be used at times and on modes for which the physical energy density 
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28| and we have tried to present 



is transplanckian. Opinion on these issues is divided 
both sides. It is worth pointing out that stress tensor fluctuations from very early times 
would also induce significant non-Gaussianities if one were to compute them perturbatively, 
using known physical laws, as we have done for the 2-point correlator. 
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Appendix A: Some relations involving spatial perturbations 

In this appendix, we will derive some relations relating to the spatial parts of the metric 
perturbations which are used in Sect. IIIIB2[ The Einstein equations, Eq. ( I40l) . may be 



26 



expressed as 






d,^d,^ + g,.V{^)\ = l/c^T^r • (Al) 



It is convenient to remove the scale factors from the conformal stress tensor by defining 

U = Tr' = fu , Tr' = afu , T^p' = a%. (A2) 

Note that T^^, are the components of the conformal stress tensor in a local orthonormal 
frame defined by dt = dt and dx'^ = adx\ The ha and hij defined in Eq. (|35|1 are the metric 
perturbations in this frame. 

Because the conformal stress tensor is itself a first order perturbation we can express its 
tracelessness using only the zeroth order metric, 

^p<xy^conf^Q ^ f,, = fu + 0{h''). (A3) 

Similar simplifications can be made to the relations implied by stress-energy conservation, 

^dtia^fu) = -ftk,k-Hfkk + 0{h^), (A4) 

a-^ a 



^dtia'^fu) = f,k,k + 0{h^). (A5) 

The expansion of the time-time component of Eq. ( jAip was performed in Sect. IIIIAI so 
here we focus on the remaining components. The first order expansions of the required 
components of the Ricci tensor are, using the metric of Eq. fl35|) . 

Rti = hk[i,k] + -ht[k,i]k + i^H^ + H)ahu + 0{h^) , (A6) 

+hk{i,j)k — -^hij^kk — -^h^ij — HaSijhtk,k + TzHa Sijh + 0{h ) . (A7) 



The remaining first-order Einstein equations become 



hk[i,k] + 2Hdi^ = -K^afu , (A8) 
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and 



d'i + mdt - \ 


h 


1 r 


-25 ij 





i^ikfkj ' i^jkfki 



l^,ij 



+ HSijh 



Bij 



2H 



6,,V' + d,dj 



B = k'T. 



V 



(A9) 



To understand what Eqs. (lASP and flA9p imply, it is useful to make a decomposition of 



hij into irreducible representations of the rotation group, 



K, ^ h^j + hi + hi - 1 (5., -3^y+l (5. 



didi\,r 1/. ^i^j\j^ 



(AlO) 



This is a decomposition into transverse-tracefree (TT), transverse (T), longitudinal (L) and 
trace parts. Here hj^ = = hfjj and hf,^ = as usual. We can make similar decompositions 
of the conformal stress tensor. 









(All) 
(A12) 



The following identities facilitate extraction of the longitudinal and trace parts. 



didj 



2 V '^ V2 



*"^J\V72 



V^ 



(A13) 
(A14) 



Equivalently, the longitudinal part is obtained by the action of the projection operator 

di d. 



'ij Y72 ' 



(A15) 



SO that h^ = Lij h^j. The longitudinal part of (IA9P is, 
{d^ + 3Hdt + \)h''+ {Hdt ~\)h 



-{2H+12HH]B + ^[2dt + 4:H]B = K^T"^ . 



(A16) 



Note that stress-energy conservation ( JASp implies, 

{dt + 3H) {aTl") = T^ 



(A17) 
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and Eq. (lASP implies 



h-AHB + n'aT, 



2„rrL 



(A18) 



^h'^+idl + AHdt-^^h 

-AHB - (qH + 2AHH)b + ^[2dt + 4H)b = . 



(A19) 



Now eliminate the h terms using Eq. ( 157|) and eliminate the resulting d^B term using 
Eq. (!58|) . The resulting simplification of (JAIOP is, 

-^(h^ -h + AHB) + (-8H^ - 127/2 ^ 4/^^ ^ _ 4^^ = /t^ f/ . (A20) 

At this point, it is convenient to switch to the variables ijj and E defined in Eq. fl9T|) . in 
terms of which 



h^ = -2^-2^-E. 



(A21) 



and 



Now Eq. (]A20p may be written as 

^2 



ip-HB + 



h = -Qtp-2^-E. 



.V^o 



K 



iJ5 + 2H— + 3H^-H] B + —U 



0. 



The time derivative of Eq. (IA22P is 



(A22) 



(A23) 



h = -6^ - 2^{E - 2HE) . 



(A24) 



Next substitute this relation and the time derivative of Eq. (1A23P into Eq. (lFr|) . and eliminate 
the dfB term using Eq. (158|) . The result is 



^{E - 2HE) =B+(^ + 3h\b- k^^ (sHU + ^U 



(A25) 
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